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ABSTRACT 

We calculate the one-point probability distribution function (PDF) for cosmic 
density 6 in non-linear regime of the gravitational evolution. Under the local 
approximation that the evolution of cosmic fluid fields can be characterized by the 
Lagrangian local dynamics with finite degrees of freedom, the analytic expressions 
of PDF are derived taking account of the smoothing effect. The validity and 
the usefulness of the local approximation are then discussed comparing those 
results with N-body simulations in a Gaussian initial condition. Adopting the 
ellipsoidal collapse model (ECM) and the spherical collapse model (SCM) as 
Lagrangian local dynamics, we found that the PDFs from the local approximation 
excellently match the simulation results in the case of the cold dark matter initial 
spectrum. As for the scale-free initial spectra given by P{k) oc k"', N-body 
result suffers from spurious numerical effects, which prevent us to give a detailed 
comparison. Nevertheless, at the quality of N-body data, the model predictions 
based on the ECM and the SCM quantitatively agree with N-body results in 
cases with spectral index n < 0. For the index n > 0, choice of the Lagrangian 
local dynamics becomes crucial for an accurate prediction and a more delicate 
modeling is required, however, we find that the model prediction based on the 
ECM provides a better approximation to the N-body results of cumulants and 
PDFs. 

Subject headings: cosmology: theory — dark matter — galaxies:clusters:general 
— large scale structure of universe — methods: analytical 

^Department of Physics, School of Science, University of Tokyo, Tokyo 113, Japan. 

^Research Center for the Early Universe(RESCEU), School of Science, University of Tokyo, Tokyo 113, 
Japan. 



- 2 - 



1. Introduction 

The probability distribution function (PDF) of the cosmological density fluctuation is a 
fundamental statistical quantity characterizing the large-scale structure of the universe. In 
a standard picture of cosmic structure formation based on the cold dark matter scenario, 
the gravitational evolution of the dark matter distribution plays an essential role for the 
hierarchical nature of observed luminous distributions. Usually, the evolution of dark matter 
distribution is believed to be developed from a small initial fluctuation with Gaussian random 
distribution. While the PDF of density fiiictTiation retains Gaussian shape in a linear regime, 
the deviation from Gaussian distribution becomes significant in the non-linear regime of 
gravitational evolution. 

A number of studies in quantifying the non-Gaussian properties of density field have 
been developed theoretically and observationally. From the numerical and the observa- 
tional study, a systematic analysis using the cosmological N-body simulation or the observed 
galaxy distribution yield various phenomenological prescription for the density PDF in the 
non-linear regime (e.g., Saslaw & Hamilton 1984; Hamilton 1985; Gaztahaga & Yokoyama 
1993; Ueda & Yokoyama 1996). Among them, the lognormal distribution has been long 
known to fit to the simulations quite accurately (e.g.. Coles & Jones 1991; Coles, Melott, 
& Shandarin 1993; Bernardeau & Kofman 1995; Taylor & Watts 2000). Recently, Kayo, 
Taruya & Suto (2001) critically examined this issue using the high resolution N-body simu- 
lation with Gaussian initial conditions, and found that the accuracy of the lognormal model 
remains valid, irrespective of the nature of initial spectra. The weak dependence of the initial 
spectra was later investigated using the phenomenological models with dark halo approach 
(Taruya, Hamana & Kayo 2003). 

On the other hand, from the analytical study, a perturbative construction of the PDFs 
has been exploited by Bernardeau (1992. 1994a) employing a field-theoretical approach and 
the predictions including the smoothing efi'ect excellently match the N-body simulations in 
the weakly non-linear regime. Beyond the perturbative prediction, however, no exact treat- 
ment is present and the non-perturbative approximation or the phenomenological approach 
taking account of the empirical simulation results are necessary. Fosalba & Gaztafiaga 
(1998a) and Scherrer & Gaztanaga (2001) proposed to use a spherical collapse model as a 
non-perturbative approximation to predict the higher-order moments and PDFs. In their 
treatment, one assumes that the Lagrangian dynamics of the local density field is sim- 
ply described by spherical collapse model. Although this approximation clearly misses the 
non-locality of the gravity in the sense that the evolution of the local density field can be 
determined by the one-to-one local mapping, the advantage of this treatment is that one 
can easily calculate the higher-order correction of the moments and PDFs. Further, it turns 
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out that the spherical collapse approximation exactly recovers the leading order results of 
perturbation theory. 

Recently, we generalize the idea of spherical collapse approximation to the local approx- 
imation in which the evolution of local density field is characterized by the Lagrangian local 
dynamics with finite degrees of freedom (Ohta, Kayo & Taruya 2003). As a demonstration, 
the PDFs were computed using the ellipsoidal collapse model. In the ellipsoidal collapse 
model, the local density at a position is expressed as the multivariate function of initial 
parameters, i.e., the principal axes of the ellipsoid given at the same position. Thus, the 
relation between the initial and evolved density field cannot be described by the one-to-one 
local mapping. As a consequence, the local approximation with ellipsoidal collapse model 
successfully explains the stochastic nature seen in the simulation, i.e., the joint probability 
between the initial and the evolved density fields, as has been reported by Kayo, Taruya 
& Suto (2001). In addition, the leading-order results from the eUipsoidal collapse model 
correctly reproduce those obtained from the exact perturbation theory. 

In the present paper, we extend the previous study to the quantitative comparison be- 
tween the local approximation and the N-body simulations. Evaluating the PDF of local 
density fields taking account of the smoothing effect, we consider the validity and the limi- 
tation of the local approximation with spherical and ellipsoidal collapse models. The PDFs 
from the spherical collapse model were previously compared with the N-body simulations in 
the case with cold dark matter (CDM) power spectrum (Schcrrcr & Gaztafiaga 2001). In 
this paper, taking account of the smoothing cfi^cct, the N-body results with the scale-free 
initial spectra as well as the CDM spectrum arc compared. 

This paper is organized as follows. In section 2, we start to review the local approxi- 
mation of one-point statistics developed by Ohta, Kayo & Taruya (2003) and briefly show 
how to compute the PDF and the moments from the Lagrangian local collapse model. As 
representative models of Lagrangian local dynamics, the spherical and the ellipsoidal col- 
lapse model are considered. Then, we consider the smoothing effect and discuss how to 
incorporate it into the model predictions. Based on this, the perturbative calculation of 
cumulants up to the two-loop order is presented and the qualitative behaviors of the model 
prediction arc discussed in section 3. In section 4, the validity and the usefulness of the local 
approximation for one-point statistics is investigated by comparing the PDFs and cumulants 
from the local collapse models with those obtained from the N-body simulations. Finally, 
section 5 is devoted to discussion and conclusions. 
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2. One-point statistics from the local collapse model 

In many analytical works on the gravitational evolution of density distributions, the cold 
dark matter distribution is often treated as the pressureless and non-relativistic fluid. This 
treatment is not exact, but in a statistical sense, it would provide a better approximation 
if the scale of our interest is large enough, where no shell-crossing appears in the smoothed 
density fields. Denoting the mass density and velocity field of fluid by p and v, the evolution 
equations for the fluid in a homogeneous and isotropic background universe are expressed as 
follows: 

f + -V-[(1 + (5H=0, (1) 
at a 

dv I 1 

— + Hv+-{v-V)v^ — V(/), (2) 
at a a 

where 5 is the density fluctuation, 5 = {p — Pm)/Pm- The quantity a is the scale factor of 

the universe and H is the Hubble parameter given by H = d/a. Gravitational potential is 

determined by the Poisson equation: 

VV = '^7^Gpma^S. (3) 

Below, owing to the local approximation, we consider the one-point PDF of the density 
fluctuation 5 and the higher-order moments taking account of the smoothing effect. 



2.1. Local approximation 

As mentioned in section 1, the analytical treatment of the one-point PDF P{S) governed 
by the fluid equations (l)-(3) is generally intractable due to the non-linearity and the non- 
locality of the gravity. Beyond the perturbative prediction, a non-perturbative treatment 
should be exploited. The local approximation is one of the way to treat the one-point PDFs 
analytically by reducing the Lagrangian dynamics of the fluid motion given by (l)-(3) to 
the local dynamics with flnite degrees of freedom. In this treatment, the time evolution of 
the local density fleld 5 at a given position is determined by the local dynamics with the 
initial condition given at the same position in Lagrangian coordinate. Thus, the solution 
of local density field can be obtained by solving a couple of ordinary differential equation 
and expressed as a function of initial parameters A = (Ai, A2, ■ ■ ■ , A„) given at a Lagrangian 
position and time, i.e., 6 = f{X,t). In this paper, we consider the spherical and the ellipsoidal 
collapse models as representative example of the Lagrangian local dynamics (see section 2.1.1 
and 2.1.2). In this case, the initial parameters of Lagrangian local dynamics correspond to 
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the linearly extrapolated density fluctuation 5i for the spherical collapse model and the 
principal axes of initial homogeneous ellipsoid, (Ai, A2, A3) for the ellipsoidal collapse model. 

Once provided the functional form of the local density /(A, t), the one-point PDF of 
the local density fleld P{6; t) (in Eulerian space) can be analytically obtained. With a slight 
modification of the definition of density fluctuation S so as to satisfy the normalization 
condition and the zero- mean of S (Ohta, Kayo & Taruya 2003), one has 

1 /■ " 

^('^'^) ^ TT^y U'^^PiW^dIS - g{\t)], (4) 

where the function g is given by 

S = g{\,t)^NE[l + f{\,t)]-l ; NE{t)^ J f^dX^^^^j^. (5) 



In equations (4) and (5), the function Pr(A) is the probability distribution of the initial 
parameters, which characterizes the randomness of the mass distribution. Prom (4), one also 
calculates the moments of the density fields: 

f r 
(5^> = / 5''P{5,t)d5 = / l[d\^PriX). (6) 

The expressions (4)- (6) are the heart of the analytical treatment in local approximation, 
which are rigorously derived by considering the evolution equations for the one-point PDPs 
(Ohta, Kayo & Taruya 2003). 



2.1.1. Spherical collapse model 

As simple Lagrangian local dynamics to calculate the function f{\,t), let us first con- 
sider the spherical collapse model (SCM). In the SCM, the evolution of local density at a 
given position in Lagrangian space is determined by the mass M inside a sphere of radius R 
collapsing via self gravity: 

-^--^ + -,R ; M=-pi?3 = const, (7) 

where A is cosmological constant. The above equation is re-expressed in terms of the local 
density defined hy 5 — (aRo/R)^ — 1: 

d'^S ^d5 4 1 fd6Y 3 
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where is the density parameter given by = SnGpm/i^H'^)- Note that the SCM can 
be regarded as the monopole approximation of the fluid equations neglecting the shear and 
vorticity (e.g. Fosalba & Gaztafiaga 1998a), since one obtains the following equation from 
(l)-(3) with a help of the Lagrangian time derivative, d/dt = d/dt + v/a ■ V: 

2 



d'^S ^dS 4 1 

h 2H 

df^ dt 31 + 5 



d5_ 
'dt 



(9) 



Uij — 



dvi ^ dvi 



2aH \dxj ' dxi 
1 / dvi dvj 
2aH \ dxi dxi 



(10) 
(11) 



where (Tij and uiij respectively denote shear and vorticity tensor. 



The exact solution of the equation (7) is obtained in the case of the Einstein-de Sitter 
universe {Vt^ = 1, A = 0) and is expressed as the function of linearly extrapolated density 
fluctuation 5i in parametric form: 

2/3 



for 5; > 0, and 



9 (?7 — sin 77) 
2 (1 — cos 77)3 



9 (sinhr/ — 77)^ 
— —- 



- 1, 



1, 



5i 



- (77 — sin 77) 



-(sinh 77 - 77) 



2/3 



(12) 



(13) 



2 (cosh 7/ — 1)3 

for 5i <Q. The relation between 5 and 5i in the above equation is fairly accurate even in the 
non Einstein-de Sitter universe (e.g. Nakamura & Suto 1995; Fosalba & Gaztanaga 1998b). 
We will extensively use (12) and (13) for later analysis. Note that in computing the density 
PDF, the linearly extrapolated density 5i should be regarded as initial parameter and is 
treated as a random variable. Assuming the Gaussian initial condition, we have 

PA^i) = -^e-(^'/'^')^/^ (14) 



27r(7; 



where the variable ui is the rms fluctuation of hnear density fleld, ui — {SfY^'^. 



2.1.2. Ellipsoidal collapse model 

The eUipsoidal collapse model (ECM) is an extension of the SCM taking account of the 
non-sphericity. In this model, the dynamics of the local density is described by the self- 
gravitating uniform density ellipsoid characterized by the half length of the principal axes 
ai {i = 1,2,3). 
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The local density of ECM is given by 

5^— 1. (15) 

According to Bond & Myers (1996), the evolution equations of half length of axes become 

—ai = -ai - 47rG ai I — ^ + K^t,i I , (16) 

r dr 2 

bi = aia2a3 / — — --^ - -. (17) 

The quantity Aext,i mimics the effect of external tidal shear, which was introduced for the 
consistency with Zel'dovich approximation: 

{Aj ^— — — - ; linear external tide , 
5 (IS) 
- bi ; non-linear external tide , 

where A^ denotes the initial perturbation of principal axis and evolves as Xi{t) — D{t)Xi{to) 
with the variable D being the linear growth rate. Note that the inclusion of the external tidal 
term is necessary to reproduce the Zel'dovich approximation when we linearize the evolution 
equations.^ In this paper, both cases of the external tidal term are considered in order to 
reveal the model dependence of the prediction, but in comparing with N-body simulation, the 
ECM results with linear external tide is presented for brevity. In presence of the external tide, 
the initial condition is specified by the Zel'dovich approximation. Identifying the variables 
Aj with the initial parameters of a^, we have 

ai{to) = a{to)[l - Xi{to)] (19) 
j^aM = d{to)[l-Xi{to)]-a{to)Xi{to). (20) 

Note that the initial parameters Aj are related to the linearly extrapolated density fluctuation 
6i by 5; = Ai + A2 + A3. Hence, the parameters Aj should be treated as the random variables. 
For the Gaussian initial condition, the distribution function of Aj is analytically expressed 
as follows (Doroshkevich 1970): 

P'i^^) = ^^^^P + 15^) (^1 - ^2)(A2 - A3)(Ai - A3), (21) 



^If this term is dropped, a consistent calculation with initial distribution (21) is also impossible. 
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where we define h — Xi + X2 + A3, I2 — A1A2 + A2A3 + A3A1. 

Note that in contrast to the SCM, the ECM can be regarded as an approximation of 
the fluid equations taking account of the effect of tidal shear but neglecting the vorticity 
(Ohta, Kayo & Taruya 2003). The equation (16) with (17) is re-expressed in terms of the 
local density S (c.f., eq.[9]): 

1 / dj di do d('^\ r 
fTi,- = — 3- - — - — - — 5ij. 23 
3H V a, «! ^2 «3 / 



2.2. Smoothing effect 

When we evaluate the statistical quantities from the N-body simulations, the smoothing 
procedure is often employed to remedy the discreteness of the particle data. In this sense, 
the smoothing effect is crucial and should be incorporated into the theoretical prediction. In 
this paper, we adopt the top-hat smoothing and the smoothed density PDFs are computed 
from both the local approximation and N-body simulation: 

S{x;R) ^ jd^y6{y)WTH{\y-x\;R); (24) 

, r 3/4rf r<R 

= jo otherwise' ^^5) 

where Wth('"; R) is the top-hat smoothing kernel of the radius R. 

A systematic method to compute the analytic PDFs taking account of the smoothing 
effect was first considered by Bcrnardcau (1994a) based on the perturbation theory. Later, 
his method was extended to the non-perturbative calculation of the PDF using SCM (Fosalba 
& Gaztanaga 1998a). We briefly review the method by Fosalba & Gaztafiaga (1998a). 

First notice the fact that in the case of the top-hat smoothing, the leading-order results 
of cumulants for local density field is not affected by the smoothing in Lagrangian space 
(Bernardeau 1994a). Extrapolating this result to the non-perturbative approximation with 
SCM, one can approximate the evolved local density with top-hat smoothing by 

S^f{Si,L), (26) 

where 5 is the smoothed density and ^^^l is the linear density fluctuation in Lagrangian space. 
To relate the quantity with that in Eulerian space, 5i^e, we recall the fact that the radius 
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R defined in the Eulerian space roughly corresponds to the radius -R(l + 5)^/^ in Lagrangian 
space. Thus, we have 

MR) (ri[R{i + Sy/^]' 

Substituting (27) into (26) and identifying the linear fluctuation Si^e with the smoothed 
linear fluctuation 6i, the relation between 5 and 6i becomes 



(2^ 



The above relation can be further simplified by the running index 7 defined by 7 = 0? log af/ d log R: 



f{5i) = f{ l + f{5i) 



7/6 



a log it 



(29) 



where n is the index of the initial power spectrum, P{k) oc A;". Notice that the above 
expression is still valid in the non power-law cases of the initial power spectrum. The 
remarkable feature in the relation (28) or (29) is that the leading order result of the cumulants 
obtained from the perturbation theory is exactly recovered by the local approximation with 
SCM. 

Owing to this noticeable fact, we extend to use the result (29) to the local approximation 
with the ECM. In this case, the smoothed density field is related to the top-hat filtered 
principal axis Aj given by: 



Adopting this form, the PDF (4) becomes 



/ l + /(A,t) A,t , 



7/6 



(30) 



1 P ^ 

P{5;R) \[dKP,{\)5D [5-g{\t)\ , 



(31) 



with ^(A, t) being 



g{\t)^NE[l + f{X),t]-l ■ Mt) ^ J l[dX,- 



(32) 



The above equations apparently seem difficult to evaluate because of the implicit relation of 
the functions / and A. However, this apparent difficulty can be eliminated by introducing 
the following variables: 



l + f{X,t) 



7/6 . 



A,;. 



(33) 
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Then the equation (31) becomes 



3 

P{6; it:) = / n dKPi {[1 + /(A', t)r^' A'} Sn [s - g{X', t) 



(34) 



where the quantities g and Ne can be recast as 



^(A',t)=iV^[l + /(A',t)]-l ; NE{t) 



l[dX 



P,{[l + /(A',t)]-^/^A'} 



l + f{X',t) 



dX' 



(35) 



Also, the Jacobian \d\j/d\'i.\ is calculated using the relation Aj = [1 + /(A, t)] '^^^X'^ as: 

3 



dX' 



[l + /(A',t)] 



-7/2 



1 - 



7 1 ^.,df 



(36) 



Note that the term —1/(1 + f )Y^i=iK{^fl'^'^'i) related to the velocity divergence and 
is expressed as 6' = '^iOii/ {Ha,i) — 3 in the case of the Einstein-de Sitter universe. Thus, 
provided the solution of the equations for ECM (eqs.[15]-[17]), the smoothed density PDF 
can be numerically evaluated by substituting the solution f{\,t) into the above equations. 
Similarly, the higher-order moments of the local density becomes 



<5^P(5; R)dS = / n dKy-^P. [(1 + f)''^'X'] 



dX[ 



(37) 



In what follows, we simply denote the smoothed density field 6 by 6. 



3. Perturbation theory in Ellipsoidal Collapse Model 

Before proceeding to compare the theoretical models with N-body simulations, it is 
useful to examine the perturbativc analysis of the local collapse models. In this section, based 
on the ellipsoidal collapse model, we present the perturbativc calculation of the cumulants 
of density field. The differences between the model predictions as well as the qualitative 
behaviors are also discussed in section 3.3. 

Here and in what follows, we treat the evolution of local density in an Einstein-de Sitter 
universe. In this case, the linear growth rate D is simply proportional to the scale factor a. 
Let us denote the half-length of the principal axis by 



ai{t) = a{t){l - CMt))}. 



(38) 
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Then, regarding the scale factor a as time variable, the equation (16) is transformed to 
In terms of the variable d, the quantities bi and S are expressed as 

/•oo J n 

^ (^-c^)(^-c^)(^-^3)/ (a-,)-^.}n.((i-o)-+r}v- -3 

= ^{30-(Ci + C2 + C3)} + fe^, (41) 
5 = C1 + C2 + C3 + A. (42) 

Note that both of the quantities bi and A are of the order of 0(C^)- 

Below, we separately give the perturbation results in the eUipsoidal collapse model with 
non-linear external tide and with linear external tide. Note that the leading-order results 
for the cumulants in each model exactly coincide with each other and reproduce those of the 
exact perturbation theory (e.g., Bernardeau 1994a,b). 



3.1. Ellipsoidal collapse model with non-linear external tide 



In the case of the model with non-linear external tide, the quantity Aext,i is given by 
(5/4)6j and the right hand side of the equation (39) becomes 3Ci/2 -|- C(C^)- Then We have 



2^2^ 3 dQ 3 

a —-^ H — a— d 

da^ 2 da 2 



i(A - 5Ci) + Ci)biS + f - Qb. 



(43) 



This equation can be perturbatively solved by substituting the series expansion Q — X]j=o ^i^^^'^ 
into the above equation. Under the initial condition (19), one formally obtains the expression 
of the coefficient ^■"'^ as 



d^ 



(2j + 3)(j-l)j! daJ 



A-6Ci + ^{l- QbiS + f[bi- Cih) 



(44) 

; J > 1, (45) 



a=0 



where we have only considered the growing mode of the solutions. Note that the coefficient 
^i^^ is the variable of the order of C(A-^). 

Based on the result (45), the pcrturbative expansion for the evolved density, 6 = 
can be constructed from the relations (15), (38) and (40). The results up to the seventh 
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order become 



<S<" = S,, (46) 

= ^^'+^^" (^^) 
,3, ^ 341 338 ^ 

567 ^ 945 3969 ' ^ ' 

„ ^ j580^,.^ig52g^,,v,^ j3408^,,j^^ 429728 

130977 ' 1091475 ' 4584195 10696455 ^ ^ ^ 

213662 . 292398464 . ^ 64182728 ,0 . 

729729 ' 638512875 ' 893918025 ' 
6541246 2 828974992 

+ 59594535'^"^^ + 96364363095'^''^'' ^^"^ 
_ 21129781 g 15739030628 4 38380501904 3 

~ 107270163 ' ^ 37246584375 ' ^ ^ 469306963125 ^ ^ 



334168450808 2 .2 250313183728 
1825082634375 ' ^ ^ 9368757523125 ' ^ ^ 



63778345006048 ^3 2272657750768 ^2 

7673012411439375 ^ 4603807446863625 ^' ^ ' 

83411812 ,7 42267062029204 ,. 63224677073056 ,4, 

5^^> = Sj H dfJi H 6fJ2 

639441621 ' 116564878828125 ' 769328200265625 ' 

3019799334120902 . 9581060236980928 , , 

H SfJf H 5, Ji J2 

12565693937671875 ^ ^ 193511686640146875 ' 

13515215809239748 , . 276922264619192 , . 

451527268827009375 ' ^ 162549816777723375 ^ ^ 
9155185965341512 . r 

^3521912696850673125 ^ ^' 

where 5/ denotes the linear fluctuation given by 5; = Ai + A2 + A3. Here, we introduced the 

quantities Ji = + xy + and J2 = (x — |/) (2a; + y){x + 2y) with the variables x and y 

being x = Ai — A2 and y = X2 — A3, respectively. 

Once provided the perturbative solution for the non-smoothed density field, cumulants 
for the smoothed density is calculated as follows. From the perturbative inversion of the 
relation (30), the smoothed density 6 is obtained and the normahzation factor Ne is first 
calculated by substituting this into the definition (32). Using the probability distribution of 
initial parameter (21), the resultant expression becomes 

, 1 2 f 10844 79 31 2 1 4 

Ne ^ I 7(7, + 7 Y 7 o-/ 

6 ' ' V 848925 4410 378 27 / ' 



/ 3891599696 _ 1248901 _ 47093 , 
^ V 511023137625 ~ 278107830 ~ 415800 
62341 19 , 1 



317520 ' 168 ' 48 



7^ - i^^' - 7^7^ (53) 
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up to the order 0{af). Then, the moments for the smoothed density (S^) is evaluated 
from the relation (S^) = J Yli g^/{^ + g) Pi{^) (cf. eq. [6]). Finally, the perturbative 
correction for the variance a"^ = (6^), the skewness = {S^)/a'^ and the kurtosis 5*4 = 
(((5^) — 3a^)/a^ are obtained and can be summarized as series expansion of af as follows: 



(J 

^3 
S4 



^2 I „ 4 



■52,6 0"f + 52,8 O'f + 
2 I c -4 



'S'3,0 + >S'3,2 +S'3,4(7; + 
5*4,0 + 54,2 O'f + 5*4,4 + 



The resultant expressions for the coefficients S2,i, ^3,1 and 54,j become 



439 167 



11 



S2,4 = 
S2,6 = 
^2,8 — 

for the variance. 



1 7 H 7" 

245 126 ' 36 ' 

3143785639 15856223 



-7 + 



55273 



1835 



127 



r7 



77 



695269575 ' 2037420 ' ' 10584 ' ' 1134 ' 648 
7932609222047169799 2321384486861437 1062497682871 



537532462889296875 



+ 



54752479031250 



-7 + 



+ 



16268385923 
495093060 



-7^ + 



61875775 
5143824 



13831 
5832 



20858087250 

6877 6 

34992^ 



-7 



53,0 = 
53,2 = 
53,4 = 



for the skewness and 
60712 

54,0 
54,2 
54,4 



34 

1041064 21946 



101871 2205 
161751288183332 

3041804390625 
11299781 



7 



415 2 10 3 

7 H 7 

126 ' 27 ' ' 



+ 



363349617641 



3476347875 



-7 + 



415283963 
5093550 



-7 



+ 



357210 



2975 
486 



■7^ + 



1841 
3888 



7 



62 7 2 

— 7 H — 7 
3 ' 3 ' 



1323 

941370178286 



+ 



1518808496 



-7 + 



18161033 



3935 



-7' 



3 



1549 



-7 



3476347875 4584195 ' 119070 ' 126 ' 648 
30144942925392628918 129392230965050887 397096017904379 o 

+ 7 + ^00^.0^^^^. 7 



13782883663828125 



27376239515625 



+ 



30126971437 
15002820 



1142621801 
2143260 



-7^ + 



6127195 
81648 



93861392625 
102005 



23328 



■7 



(54) 
(55) 
(56) 



(57) 
(58) 

(59) 

(60) 
(61) 

(62) 

(63) 
(64) 

(65) 



for the kurtosis. 
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3.2. Ellipsoidal collapse model with linear external tide 

For the pcrturbative solution in the model with the linear external tide, the calculation 
is slightly reduced if we introduce the quantities A = C1 + C2 + C3! -B = Ci~C2 and C = C2 — Cs- 
Recalling the fact that the external tidal term becomes Aext,i = K ~ (-^i + A2 + A3)/3, the 
evolution equation (39) is rewritten with 



" V + i-Tta - 2-^ = 2A-2^^-iEC.(* + 2A.,.), (66) 



i=l 

rfCpB 3 dB 3,, ,, 3^,,_ ^x,-, 1,-^ 

" V + 2« Ifa - 2'^' - = i^{(l-C*-(l-C*}-2*« 

3 

— 2(Cl'^ext,l — C2Aext,2), (67) 

3 

— 2(C2Aext,2 — CsAext.s)- (68) 

Note also that the right hand sides of the equation (66)-(68) are of the order of 0{\^). 

Similar to the procedure in section 3.1, the perturbative expansion for the density 5 is 
constructed from the perturbative solutions of A, B and C. After a tedious but a straight- 
forward calculation, the perturbative solutions up to the seventh order become 

5^'^ = 5u (69) 

,,3^ 341 3 1538, ^ 4 ^ 

= 567^' + 4725^'^^ + 405^- ^''^ 
_ 55805 , 952144 345088 12368 , 

^ ~ 130977^' + 2480625'^' + 16372125'^''^' + 363825"^^' ^^^^ 

213662 . 237342074 ,3 , 93363344 t 

5^^> = (5f H SfJi H 5f J2 

729729 ' 621928125 ' 3192564375 ' 

52865818 , ,0 135052 ^ ^ 

H 5, Ji H Ji J2 , (73) 

638512875 ^ 34827975 ^ ^ 

(6) _ 21129781 g 73816004896012 4 29134959410408 ,3 

107270163 ' ^ 215099024765625 ^ ^ ^ 879950555859375 ^ 
45534497984 „ 1416570594232 , ^ ^ 

355535578125 ' ^ 129059414859375 ' ^ ^ 

797014912 . 162352 „ 

H -^1 H -^^-J, (74) 

132368630625 ^ 1578740625 ^ ^ ' 

83411812 ,7 24700151148166244,. 
639441621 ' 85566392051765625 ' ^ 



49159400006961656 ^ 18569055254261594 

^1480956785511328125 ' ^ ^ 116681443706953125 ' ^ 
13101172588796 ^ r 9323215177292 , . 

684531136414125 ' 488950811724375 ^ 

1452480606736 , ,0 11003633175272 . ^ 

4400557305519375 ' ^ 10267967046211875 ^ ^' 



The normalization factor of the PDF is 



/ 69668 701 31 o 1 4 \ 4 



^ 6^*^' ^ I 3898125 ^ 198450^ 378^ 27^ j 



/ 5033872069084 3873942169 38616157 
145297074296875 ^ 2234-; 

286169 , 19 . 1 



^ V 645297074296875 ^ 223479506250^ 509355000^ 



1587600 ' 168 ' 48 

Then, the coefficients of the perturbative correction for cumulants become 

57137 167 11 2 

^^''^ ~ 33075 ^ 126^ ^ 36^ ' 

469828713881 17130160379 488945 2 1835 3 127 4 

~ 111739753125 ^ 2292097500 ^ ^ 95256 ^ ^ UzV ^ 648^'' ' 

996244294855051546571 152613969392185373 683964582869801 

~ 74870593045294921875 ^ 3871782445781250 ^ ^ 14079208893750 ^ 
1980638022487 3 306289019 , 13831 , 6877 



H 7 H 7 H 7 H 7 

61886632500 ' 25719120 ' 5832 ' 34992 ' 



for the variance, 



34 

'53,0 = y + ^' 

646404856 327062 415 2 10 . 

Sr.n — 1 1 H 1 H ^ 

63669375 33075 ' 126 ' 27 ' ' 
77881923244216108 80186055186641 2052918391 2 

^'^ ~ 1505693173359375 ^ 782178271875 25467750 ^ 

56239289 3 2975 , 1841 , 

H 7 H 7 H 7 

1786050 ' 486 ' 3888 ' 

for the skewness and 

60712 62 7 2 

'° " T323" ^ y ^ ^ 3^ ' 

210688932175742 188859083824 90600877 , 3935 3 1549 . 

Sa2 = \ 7 H 7 H 7 H 7 

782178271875 573024375 ' 595350 ' 126 ' 648 ' 

102133149992759420855618 7007053550215029257 

^'^ ~ 47644922847005859375 ^ 1505693173359375 ^ 

29469590927547677 2 13726819977457 3 

^ 7039604446875 ^ ^ 6876292500 ^ 

5700735749 . 6127195 , 102005 « 

H 7 H 7 H 7 

10716300 ' 81648 ' 23328 ' 
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for the kurtosis. 



3.3. Differences between model predictions 

To understand both the quahtative and the quantitative behaviors of the above two 
predictions, we here briefly discuss the systematic dependence of the perturbative results. 
Table 1 summarizes the numerical values of the coefficients, S2,i up to the three-loop order, 
and Ss^i and S^^i up to the two-loop order for each model with various spectral indices. And 
using these results, we plot the cumulants up to the two-loop order in figure 1, where the 
cumulants are normalized by the leading-order results af, 3^^ and 84^ (from top to bottom) 
and are depicted as function of linear variance o"/. The results from the SCM (short- dashed) 
are essentially the same results as those obtained by Fosalba & Gaztahaga (1998a). Note that 
the leading-order results for the cumulants in all model predictions rigorously coincide with 
those obtained from the exact perturbation theory (e.g., Bernardeau 1994b), irrespective of 
the choice of the external tidal term. 

Figure 1 shows that the differences between the model predictions are generally small 
for the spectral indices n < and these are expected to become negligible as decreasing 
n, approaching the non-smoothing results (n — —3) as obtained previously (Ohta, Kayo & 
Taruya 2003). For n — —2, the predictions up to the one-loop order give 

^ af + 0.6laf (86) 
^3 fti 3.86 + 3.21(7^2 (^37^) 

for the SCM (Fosalba & Gaztaiiaga 1998a) and 

a; af + 0.88af (88) 
S'3 ^ 3.86 + 3.18(7^ (89) 

for the exact perturbation theory (Scoccimarro & Frieman 1996; Scoccimarro 1997). Com- 
paring the above results with Table 1, the numerical values of the coefficients are close to 
those from the EGM prediction. On the other hand, for the n = case, a large discrepancy 
appears in the variance (X^ (Fig.l). The skewness and the kurtosis also show a relatively 
large difference. These behaviors are indeed consistent with the PDF shown in figure 5. The 
reason why the discrepancy in the model predictions become large as increasing n is partially 
ascribed to the Jacobian in the smoothed density PDF (34). The expression of the Jacobian 
\dXj/dX',^\ in (36) contains the quantity related to the velocity divergence 9 multiplied by the 
spectral dependent factor, 7/6. In previous study of the non-smoothing case (Ohta, Kayo & 
Taruya 2003), we found that while the differences between the model predictions are almost 
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negligible for the density fields, a large difference appears in the velocity divergence d. This 
readily implies that the differences between the model predictions also become significant in 
the present case, depending on the factor 7 = — (n + 3). In other words, for a large deviation 
from 7 = 0, the model predictions sensitively depend on the choice of the Lagrangian local 
dynamics in local approximation. That is, not only the evolution of local density but also 
the evolution of velocity field should be devised to approximate the fluid dynamics precisely. 
This point is important and should be kept in mind when comparing the prediction with 
N-body simulations (see Sec.4.2). 



4. Compcirison with N-body simulations 

Wc arc now in position to discuss the validity and the usefulness of the local approxi- 
mations comparing the theoretical prediction with N-body simulations. For this purpose, we 
specifically use the N-body data for the scale- free models with initial power spectra P{k) oc /c" 
{n — —2, —1, 0) (Jing 1998), as well as the CDM model with cosmological constant (ACDM; 
Jing & Suto 1998). The gravitational force calculation is based on the P^M algorithm. All 
the models employ A'" — 256^ dark matter particles in a periodic comoving cube I/|ox) where 
the box size of the ACDM model is chosen as Lbox = 300/i~-'^Mpc. While the scale-free 
models assume an Einstein-dc Sitter universe, cosmological parameters of the ACDM model 
are set as (fimi ^Ai o's) = (0.3,0.7,0.7,1.0), where the normalization (Xg means the linear 
rms fluctuation at top-hat smoothing radius R = 8/i~^Mpc. For the scale-free models, the 
normalization of the density fluctuation is determined by setting the linear rms fluctuation 
to unity &t R — O.ILbox- Below, we first present the results in ACDM model (Sec.4.1). 
The PDFs in the scale-free models are compared in section 4.2. In comparing the prediction 
with simulation, we also present the perturbation results of local approximation obtained in 
previous section. 



4.1. ACDM model 

The validity of the local approximation using the SCM has been previously studied by 
Scherrer & Gaztaiiaga (2001) in the case of the standard CDM model and a good agreement 
with N-body simulation was found. We thus expect that the local approximation with both 
the SCM and the ECM also provides an excellent agreement with N-body simulation in the 
case of the ACDM model. 

Figure 2 shows the PDFs obtained from the N-body data for the top-hat smoothed 
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density fields {open- squares). The error bars indicate the la errors among three different 
reaUzations. The smoothing radii are chosen as i? = 16, 8, 2^~^Mpc (from top to bottom). 
In figure 2, the PDFs from the SCM and the ECM with hncar external tide are depicted as 
long-dashed and solid lines, respectively. We also calculated the PDFs from the ECM with 
non-linear external tide, but the results are almost the same as obtained from the ECM with 
linear external tide. Clearly, these models almost coincide with each other and the agreement 
with N-body results is excellently good. For comparison, we also plot the empirical model 
of the lognormal distribution {short- dashed): 



with (7ln = log(l + (T^). Here the quantity a denotes the variance of the local density 
field, which is estimated from the N-body simulation. Albeit the simplicity of the analytical 
expression (90), the lognormal PDFs also approximate the N-body results quite accurately. 
Agreement with N-body simulation still remains good even at the high-density tails of small 
radii, R = 2 and 8/i~^Mpc. Bernardeau & Kofman (1995) discuss the successful lognormal 
fit of the PDF in the CDM models based on the perturbation results. 

In order to check the accuracy of the model predictions, we quantify the cumulants of 
the density fields. In figure 3, we plot the variance cr^, the skewness = (5'^)/a'^ and the 
kurtosis S4 = ((5^) —3a'^)/a^ as a function of smoothing radius (from top to bottom). In each 
panel, the crosses with error bars indicate the results obtained from the N-body simulations, 
while the open-squares show the results from the local approximation with ECM, in which 
the moments (5^) are calculated from the full knowledge of the PDF P{S) (see eq.[37]). On 
the other hand, the short-dashed and the long-dashed hues are the perturbative calculations 
of the cumulants based on the ECM up to the one-loop order and the two-loop order, 
respectivcly(scc Sec. 3. 2). As a reference, we also plot the leading-order(tree-lcvcl) prediction 
in dotted lines. In contrast to a naive expectation from figure 2, the ECM prediction based on 
the PDF significantly deviates from the N-body results at the smaller radius, R ^ 8/i~^Mpc, 
although it roughly match the perturbation results up to the two-loop order. 

Kayo, Taruya & Suto (2001) remarked that the origin of this discrepancy might be due 
to the fact that the density field 6 in N-body simulation does not extend the entire range 
between —1 and +00, but rather is limited in the range, ^min < d < 5max owing to the finite 
size of the simulation box. Indeed, a closer look at figure 2 reveals that there exist the sharp 
cutoff at the high-density tails of simulated PDF {arrows in each panel). To examine the 
significance of this effect, the ECM predictions of the cumulants taking account of the finite 
range [^min, (^max] are plotted in figure 3 (cutoff-1, open-triangles). The cutoff values ^min and 
(^max are estimated from the N-body data. Further, in figure 3, we plot the prediction taking 
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account of the cutoff values determined from the simple assumptions that (i) the major effect 
for the finite range of 5 comes from the finite sampling effect and (ii) the theoretical PDF 
is correct if the box size of the simulation becomes infinite. We then have (Kayo, Taruya & 
Suto 2001): 



where P{()) is the PDF obtained from the ECM. The resultant cutoff values (5min, ^max) are 
plotted as function of smoothing radii in figure 4. 

The resultant amplitudes of the cumulants are significantly reduced and the model 
prediction turns out to reproduce the simulation data very well in the case using the N- 
body data (cutoff- 1). This readily implies that the apparent discrepancy in the cumulants 
mainly comes from the hmited range of the density in PDF and one concludes that the local 
approximation with ECM provides an accurate prediction for both the density PDF and 
the cumulants in the ACDM model. Similarly, one expects that the accuracy of the model 
prediction still remains good for the local approximation with SCM. Note, however, that the 
predictions with (91) still exhibits the discrepancy at the smaller scales R ^ 4/i~^Mpc. This 
means that the finite sampling effect might be a major numerical effect for the limited range 
of the density, but still not enough. The discreteness effect could be an important source 
for causing the discrepancy on small scales. This point is particularly important for self- 
consistent calculation of cumulants and should be treated carefully. Keeping these remarks 
in mind, we next proceed to the scale-free models. 



N-body simulations in scale-free model with index n < —1 generally suffer from spuri- 
ous numerical effects compared to the CDM case, which significantly affects the statistical 
properties of mass distribution (e.g., Colombi, Bouchet & Hernquist 1996). This is not ex- 
ceptional in our case. Figure 5 shows the density PDFs for n = — 2 to models at different 
smoothing radii; R — 0.15, 0.1 and 0.02Lbox- For larger smoothing radii, the PDF obtained 
from the simulations (open-squares) has a sharp cutoff at high-density tails. Because of this, 
the amplitude of the cumulants is significantly reduced and the N-body simulation fails to 
even reproduce the leading-order results of perturbation theory (see Figs. 6 and 7). Never- 
theless, focusing on the moderately non-Gaussian tails of PDF, we find that the predictions 
based on the SCM and ECM shows a good agreement with N-body simulations for smaller 
spectral indices, n = —2, —1. Especially at the non-linear scale [R = 0.02Lbox), the PDF 
from the ECM seems to improve the prediction, compared to the results obtained from the 
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4.2. 



Scale-free models 
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SCM. This is even true for the spectral index n = 0, indicating that the ECM provide a 
more physical model of Lagrangian local dynamics than the SCM. Interestingly, the lognor- 
mal model also provides a good approximation to the N-body simulations, irrespective of 
the nature of the initial spectra and the smoothing scales. Apparently, this contradicts with 
the perturbation results which predict the strong spectral dependence in weakly non-linear 
regime (Bernardeau 1994b). As stated by Bernardeau (1994a) and Bernardeau & Kofman 
(1995), however, the perturbation results themselves resemble the lognormal PDF near the 
index n — —1. Further, a systematic comparison with lognormal prediction done by Kayo, 
Taruya & Suto (2001), who basically used the same N-body data as ours, shows that the 
lognormal model prediction tends to deviate from simulation for n = 0, +1 in the weakly 
non-linear regime even if the cutoff of the PDF is taken into account. Thus, no serious 
contradiction is found at least at the quality of our data. 

Now, consider the cumulants of the density fields. Figures 6 and 7 show the variance, 
the skewness and the kurtosis, which are compared with the prediction from the SCM and 
the ECM, respectively. As anticipated from figure 5, the amplitude of the cumulants from 
N-body data are significantly reduced and the prediction from local approximation without 
cutoff (open-squares) generally overpredicts the simulation results. A more serious aspect 
is that the simulation does not converge to the tree-level results of perturbation theory. 
Even the predictions up to the one-loop correction overpredict the simulation results. As 
previously remarked by Colombi, Bouchet & Hernquist (1996), the recovery of tree-level 
results is difficult for the limited size of N-body simulations and a care must be taken 
in order to correct the spurious numerical effects. Colombi, Bouchet & Hernquist (1996) 
devised to correct the finite volume effect to explore the scaling properties of the PDFs. 
Nevertheless, the general tendency between our simulations and those obtained by Colombi, 
Bouchet & Hernquist (1996) is quite similar, i.e., the non-linear correction to the skewness 
and the kurtosis seen in the N-body results is generally small and their amplitudes are nearly 
constant on scales. 

In figures 6 and 7, repeating the same procedure as in the ACDM case, we plot the 
cumulant predictions from the PDF taking account of the cutoff: open-triangles (cutoff- 1) 
and fiUed-triangles (cutoff-2). Then, the overall behaviors of the model predictions seem to be 
greatly improved. The agreement between the model prediction and the N-body simulation is 
reasonably good in both cutoff-1 and cutoff-2 cases except for the strongly non-linear regime, 
ai ^ 5 . A closer look at these figures shows that the prediction based on the ECM provides 
a better approximation than that of the SCM especially for the n — case. Note, however, 
that most of the prediction for the variance slightly overpredict the simulation except for 
the n = —2 case. The overprediction might be partially ascribed to the cutoff of the density 
[(^min, (^max], siucc the the width of the PDF quantified by the variance sensitively depends 
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on the cutoff. On tfie otfier fiand, tlie reduced amplitudes S3 and 5*4 ratlier cfiaracterize tfie 
sliape of tlie PDF, wliicli could be, in principle, less sensitive to the cutoff of the high-density 
tails than cumulants themselves if the tails of PDF converges to zero enough. Actually, the 
shape of the predicted PDF resembles that of the PDF in the simulations as in figure 5. The 
same tendency has been previously reported by Fosalba & Gaztanaga (1998a). 

Thus, at a level of the quality of N-body data, the local approximation with both the 
SCM and the ECM works reasonably well and the approximation with ECM even slightly 
improves the prediction. Of course, one must still care about the cutoff density arising 
from the spurious numerical effects. In this respect, the validity of the local approximation 
may just reach at an acceptable level. In addition, another important caveat is drawn 
from the strong model dependence of the predictions. As discussed in section 3.3, the local 
approximation itself becomes more sensitive to the choice of the Lagrangian local dynamics 
as the deviation from the spectral index n — —3 gets larger. This is clearly seen in the 
strong model dependence of the one-point PDF (Fig. 5) and the non-linear correction for 
the cumulant prediction (Fig. 1). Further recall the fact that the linear variance ai scales 
as ai oc This means that slight decrease of the smoothing radius R significantly 

increases the linear variance on non-linear scales, o"; ^ 1. Hence, one expects that the model 
prediction for n > suffers from the non-linear corrections more seriously, compared to the 
initial spectra n — —2, —1 or ACDM. Therefore, the local approximation with ECM should 
be used with caution for the index n > 0. 



5. Conclusion and discussion 

In the present paper, we critically examined the validity and the usefulness of the local 
approximation to the PDF and the cumulant predictions. Adopting the ellipsoidal and the 
spherical collapse model as representative model of the Lagrangian local dynamics, the PDFs 
and the cumulants are calculated taking account of the smoothing effect and the resultant 
predictions are compared with the N-body simulations with a Gaussian initial condition. 
Due to the cutoff of the density arising from the spurious numerical effects, the detailed 
comparison in cumulants becomes difficult, and the correction for the cutoff density should 
be self-consistently incorporated into the model prediction. At a level of the quality of the 
N-body data, however, the local approximation with both SCM and ECM successfully re- 
produces the N-body results for the PDFs and the cumulants, although a self-consistent cal- 
culation of local approximation presented in this paper (labeled by "cutofT-2"' ) is still needed 
to be improved. This is indeed the case of the ACDM model and the scale-free models with 
indices n — —2, — 1. For the scale- free model with n — 0, while the discrepancy between the 
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model prediction and the simulation result is manifest in the local approximation with SCM, 
the agreement with N-body results still remains good for the ECM prediction. The detailed 
discussion reveals that the prediction based on the local approximation sensitively depends 
on the slope of the initial spectrum and that the predictions for n > become more sensitive 
to the non-linear dynamics of the local collapse model. Thus, a more delicate modeling of the 
Lagrangian local dynamics is required for an accurate prediction. Taking this point carefully, 
we therefore conclude that the local approximation with SCM and ECM provides an excel- 
lent approximation to the N-body simulations for CDM and scale-free models with n < in 
both the linear and the non-linear regimes, ai ^ 5, while the local approximation should 
be used with caution in the n >0 cases. 

In this paper, we found that the predictions based on the ellipsoidal collapse model 
somewhat improve the approximation, however, the degree of the improvement is not so 
large as long as a CDM-hke initial spectrum (i.e., effective spectral index riefr = — 3 — 
(i log erf /(i log i? < 0) is concerned. Compared to the prediction from the spherical collapse 
model, the calculation of PDF from ellipsoidal collapse model is rather complicated and 
require a time-consuming numerical integration. It seems that the spherical collapse model 
provides a simpler prescription for the PDF in real space and is practically more useful than 
the ellipsoidal collapse model. However, if one considers the one-point statistics in redshift 
space, the situation might be changed drastically. As reported by Scherrer Sz Gaztanaga 
(2001), the local approximation with spherical collapse model only provides a good approx- 
imation to the redshift space PDFs when ai ^ 0.4. A part of this reason is ascribed to the 
fact that the model prediction cannot recover the linear perturbation result, referred to as 
the Kaiser effect (Kaiser 1987); the variance in the redshift space, cr^ is related to the one in 
the real space as: 



In contrast to the spherical collapse model, which leads to the incorrect prediction cr^ = 
(1 + fn/'^Ycrf, the Kaiser effect (92) can be correctly recovered by means of the ellipsoidal 
collapse model. The derivation of equation (92) is presented in appendix A. This fact 
is very interesting and also provides an important suggestion that the non-sphericity of the 
Lagrangian local dynamics play a crucial role in computing the one-point statistics in redshift 
space and is indeed essential for an accurate prediction. The detailed analysis of the model 
predictions in redshift space is now in progress and will be described elsewhere. 
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A. Derivation of Kaiser factor by ellipsoidal collapse model 

In this appendix, we derive the Kaiser effect (92) from the eUipsoidal collapse model. 
Assuming the distant-observer approximation, let us consider the local density located at 
r — (ri, r2, r^) in real space and choose the third axis as the line-of-sight direction. Denoting 
the corresponding coordinate in the redshift space by s = (51,52,53), the relation r and s 
becomes 

si = ri, S2 = r2, S3 = ra + V3/H, (Al) 

where V3 is the hne-of-sight component of the peculiar velocity field. Then the local density 
in redshift space, Ps can be expressed in terms of the quantities in real space as follows: 

dM dM p 



(A2) 



The peculiar velocity field at the position r is described by the motion of the homogeneous 
ellipsoid. Introducing the new coordinate along the principal axis of ellipsoid, r' — (r^, ri^, r'^), 
the peculiar velocity v is given by 

^i = ^i+[^.- < (A3) 

in the new coordinate system. Here, Vi is the bulk velocity of the ellipsoid. The new 
coordinate r' does not necessarily coincides with the original one r. Rather, it is related 
to the original coordinate through the Euler angle, i.e., r' = R^{'il))R2{9)R-i{(l))r , where Ri 
is the rotational matrix with respect to the i-axis. Using this fact, the quantity (A3) is 
transformed into the original frame and we obtain 



i=l ^ 



i3 



= [—-H]cos'^l^sin'e+[—-H]sm^^sin^e+{'^-H]cos'd. (A4) 



Oil / \0t2 J 

Here, we neglect the bulk velocity. Substituting the above equation into (A2) yields 

1 + 5. = ^^^ -. ^, (A5) 



cos ■V'sin^6'H sin^^sin cos 6* 

H \ai 0L2 0L2, 

with the quantity bg being the density fluctuation in redshift space. Note that the expression 
is exact under the distant-observer limit. In the linear perturbation, the quantity ctj in (A5) 
is replaced with a(l — Aj) (eq.[19]). We have 

5^ = 5; ( Al cos^ i\) sin^ 61 A2 sin^ i\) sin^ 61 A3 cos^ Q) , (A6) 



-25- 



where fn = dlnD/dlna ~ Q^. Hence, the linear variance in redshift space, is expressed 
as 

al = {SD =(^1 + ^fn^ {(X,+X,+X^f) + Lf^{3{Xl+Xl+Xl)+2{X,X,+X2Xs+XsXi)), (A7) 

where we have taken the averages over the angles, ip, 9 and (p. Finally, the ensemble averages 
over the variable Aj are taken with a knowledge of the distribution function (21): 

= (l + + Ift) af. (A8) 

This is exactly the Kaiser effect. Note that in the cases adopting the spherical collapse 
model, the variable Aj means Si/3. Hence, the ensemble average over 5i in equation (A7) 
immediately yields incorrect prediction, erf = (1 + fn/SYaf. 
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Table 1. Coefficients of perturbative correction for the cumulants of the density field 







non-hnear external tide 






linear external tide 




n 


-3 


-2 


-1 





-3 


-2 


-1 





S2,4 


1.79 


0.772 


0.363 


0.566 


1.73 


0.708 


0.299 


0.501 


"52,6 


4.52 


0.539 


0.0364 


0.360 


4.20 


0.442 


-0.0202 


0.165 


S2,8 


14.8 


0.294 


-7.50 X 10-^ 


0.158 


13.3 


0.199 


-0.0113 


-0.235 


^3,2 


10.2 


3.19 


0.525 


0.00379 


10.2 


3.19 


0.587 


0.130 


^3,4 


53.2 


4.20 


-0.0169 


0.0646 


51.7 


3.98 


0.0119 


0.246 


Sa,2 


271 


63.2 


6.67 


-0.0227 


269 


63.1 


7.32 


0.651 


'S'4,4 


2.19x10^ 


146 


0.614 


0.364 


2.14x10 


3 141 


1.16 


1.41 
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Fig. 1. — Differences between the model predictions in the variance, the skewness and the 
kurtosis. The perturbative results up to the two-loop order are shown in various models of 

the Lagrangian local dynamics. While the solid {long-dashed) lines represent the predictions 
based on the ECM with non-linear (linear) external tide, the short-dashed lines indicate the 
results obtained from the SCM (Fosalba & Gaztanaga 1998a). 
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Fig. 2. — PDFs of density field in ACDM model with top-hat smoothing window; R — 

2(bottom), 8{middle) and 16h~^Mpc{top) . The open-squares represent the N-body results. 
The error bars indicate the la variation among three different realizations. Values of ai in 
each panel are the linear variance at each smoothing radius. The arrows indicate the mean 
value of the cutoff density 5max- Note that the cutoff density at i? = 2/i~^Mpc reaches at 60f . 
Solid lines: the prediction based on the ECM with linear external tide. Long-dashed lines: 
the prediction obtained from the SCM. Short-dashed lines: the lognormal PDF adopting the 
variance calculated directly from the simulations. 
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Fig. 3. — Variance {top), skewness (middle) and kurtosis (bottom) of the density field in 
ACDM model as function of smoothing radius R. The crosses with error bars represent 
the results from N-body simulations. The open-squares show the prediction from the ECM 
with linear external tide based on a full knowledge of the PDF (eq.[37]). The open- and 
filled-triangles are the same prediction as open-squares, but taking account of the limited 
range of the density PDF, [(5min, f^max] (cutoff-1, cutoff-2). While the cutoff values in open- 
squares are estimated from the N-body data, djoin and 5max in filled-triangles are determined 
from equation (91) with a help of theoretical PDF. Long-dashed lines: the perturbative 
predictions of cumulants based on the ECM up to the two-loop order. Short-dashed lines: 
the perturbative predictions up to the one-loop order. Dotted fines: the leading order results 
of the perturbation theory. 
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Fig. 4. — The cutoff values of the density field (5min and (5max as function of smoothing radii. 
While the open-triangles show the cutoff density estimated from the N-body simulations, the 
filled-triangles represent the values obtained from the finite sampling effect (91). 
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Fig. 5. — Same as Fig. 2, but in the scale-free models {n — —2, —1,0); R — 0.02 (bottom), 
0.05 (middle) and 0.15Lbox (top). 
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Fig. 6. — The variance, the skewness and the kurtosis of the N-body results in scale-free 
models, compared with the SCM predictions. The crosses with error bars and solid line 
represent the N-body simulations. The squares show the predictions from the SCM based on 
a full knowledge of the PDF (eq.[37]). The open-, filled-triangles are the same as squares, but 
we take account of the cutoff of the density field (cutoff-1, cutoff-2). While the long-dashed 
lines represent the perturbative predictions up to the two-loop order, the short-dashed lines 
indicate the results up to the one-loop order. As a reference, we also plot the leading-order 
results of perturbation theory in dotted line. 
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Fig. 7. — Same as figure 6, but the theoretical predictions based on the ECM with hnear 
external tide approximation are presented. 



